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ABSTRACT 

Point lens microlensing events with impact parameter close to the source stellar 
radius allow the observer to study the surface brightness profile of the lensed source. 
We have examined the effects of photospheric star spots on multicolour microlensing 
lightcurves and investigated the detectability of such spots in different wavebands as 
a function of spot temperature, position, radius and lens trajectories. We include the 
effects of limb darkening and spot projection as a function of position on the stellar 
disk. In particular we apply the updated, state-of-the-art 'NextGen' stellar atmosphere 
models of Hauschildt et al. which predict very strong limb darkening, and which are 
likely to be applicable to the source stars considered here. 

Our results indicate that star spots generally give a clear signature only for transit 
events. Moreover, this signature is strongly suppressed by limb darkening for spots 
close to the limb, although the spots may still be clearly detected for favourable lens 
trajectories. 

It is also clear that intensive temporal sampling thoughout the duration of the 
transit is necessary in order for such events to be effective as a tool for imaging 
stellar photospheres. Nonetheless, with sufficiently well sampled light curves of good 
photometric precision, microlensing can indeed place useful constraints on the presence 
or otherwise of photospheric starspots. 

Key words: gravitational lensing: stars - stars: spots - stars: atmospheres - galaxies: 
stellar content 



1 INTRODUCTION 

Recently gravitational microlensing has been developed as a 
tool to probe the distribution and nature of dark matter in 
the Milky Way (see Paczynski 1996 for a detailed review), 
with several hundred candidate events detected to date. 

In the past few years, several authors have established 
that analysis of microlensing events can, in principle, pro- 
vide much useful information about the star being lensed, 
in the particular case where the star has significant angular 
extent. Such 'finite source effects' were considered initially 
for a uniform circular disk (see e.g. Bontz 1979; Nemiroff 
and Wickramasinghe 1994; Gould 1994, 1995; Witt and Mao 
1994; Peng 1997 ). These calculations demonstrated that the 
finite extent of the source would produce a deviation from 
the standard lightcurve for a point source, allowing estima- 
tion of the source radius and hence the Einstein radius of the 
lens. Heyrovsky and Loeb (1997) extended this treatment to 
the microlensing of a uniform elliptical source. 

The microlensing signature of a non-uniform disk (due 
to, for example, limb darkening) was considered by Valls- 



Gabaud (1994), Witt (1995), Bogdanov and Cherepashchuk 
(1995, 1996) and Simmons, Newsam and Willis (1995) and 
Simmons, Willis and Newsam (1995), who showed that the 
microlensed lightcurves would display a chromatic signature 
as the lens effectively sees a star of different radius at dif- 
ferent wavelengths. Further studies of the chromatic signa- 
ture of extended sources followed in, e.g. Gould and Welch 
(1996), Valls-Gabaud (1996), Sasselov (1997) and Valls- 
Gabaud (1998). These authors discussed the possibility of 
using microlensing to constrain stellar atmosphere models 
by determining the radial surface brightness profile of the 
lensed star as a function of wavelength. Hendry et al. (1998) 
investigated a non-parametric approach to inverting the ra- 
dial surface brightness profile from multi-colour lightcurves, 
using the Backus- Gilbert method (see also Gray and Cole- 
man 2002). 

The calculations of Simmons, Newsam and Willis (1995) 
were carried out for a grey model atmosphere with a linear 
limb darkening law; this simple model was improved in Valls- 
Gabaud (1998) to incorporate linear, quadratic and loga- 
rithmic limb darkening laws, with coefficients calculated for 
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the Johnson wavebands U to K. More recently Heyrovsky, 
Sasselov and Loeb (2000) have used one particular model 
atmosphere for a red giant star to explore the efficacy of 
high-precision photometric and spectroscopic microlensing 
data as a detailed probe of red giant atmosphere models. 

The study of stellar atmospheres is one of several spe- 
cific astrophysical applications of extended source events 
which have been discussed in recent literature. Others in- 
clude the measurement of stellar rotation (e.g. Maoz and 
Gould f994), diagnosis of stellar winds (fgnace and Hendry 
f999), and detection of extra-solar planets (Gaudi and 
Gould 1997). 

A subject which has received comparatively little atten- 
tion to date is the microlensing signatures of non-radial sur- 
face brightness profiles on the source disk - as would result 
from e.g. the presence of photospheric starspots. Sasselov 
(1997) presented light curves for a simple model of a cool 
circular spot on the photosphere of a red giant star. This 
treatment was extended in Heyrovsky and Sasselov (2000), 
hereafter HS00, to consider the signatures of hot and cool 
circular spots as a function of spot position, lens impact pa- 
rameter, stellar radius and spot area. The authors presented 
maps of "spot detectability" , adopting a change in the mi- 
crolensed flux of 2% as the criterion for detectable features 
in that particular model atmosphere of a red giant. Han 
ct al. (2000) extended this analysis for the case of caustic 
crossings produced by a binary lens, but did not consider 
the effects of limb darkening nor used model atmospheres. 

The stellar atmosphere model adopted in HS00 - based 
on those developed in Heyrovsky, Sasselov and Loeb (2000) - 
included the effects of limb darkening, which - as one would 
expect - suppresses the signature of a hot or cool feature 
close to the stellar limb. The numerical procedure adopted 
both by HS00 and Han et al. (2000), however, ignored geo- 
metrical foreshortening which would result in a circular spot 
appearing as an ellipse in projection, as the spot is displaced 
from the centre of the disk. 

In Bryce, Hendry and Valls-Gabaud (2002) we improved 
the geometrical treatment of HS00 and Han et al. (2000) by 
introducing an accurate model of the effects of foreshorten- 
ing for a circular spot plus an exploration of the effects of 
limb darkening within the spots themselves as opposed to 
the uniform spots in these two papers. Here we extend the 
analysis of that paper to include a more comprehensive ex- 
ploration of the effects of non-radial surface brightness vari- 
ations as a function of source, spot and lens parameters. In 
particular we extend our model atmospheres calculations to 
incorporate the recent 'Next Generation' (NextGen) atmo- 
sphere models of Hauschildt et al. (1999a, b), which include 
a more comprehensive treatment of spherical geometry and 
the effects of molecular opacity in the outer atmospheres of 
cool giants. 

The gravitational imaging of stellar photospheres is a 
potentially very useful tool for stellar astrophysics. The at- 
mospheres of cool giants are relatively poorly understood 
- not least because of their highly evolved and thus theo- 
retically complex nature. Fundamental issues such as their 
typical rotation speeds remain uncertain, and in fact spec- 
troscopic observations of microlensing events have already 
been proposed as a useful diagnostic of red giant rotation 
(Gould 1997). 

Detecting photospheric spots in ways different from, 



e.g., rotationally induced photometric modulation or 
Dopplcr imaging could bring considerable insight on their 
properties. For instance it is not clear what the covering fac- 
tor could be, particularly if many spots are present which 
bias the determination of unspotted surface properties. In 
Doppler imaging analyses the temperature of the spots is 
usually cooler than the photosphere of the stars under study, 
typically by 600 - 1200 K, but this may be at least in part 
due to selection effects. There is an interplay between spot 
area and temperature which microlensing could perhaps re- 
move. Although Doppler imaging has revealed irregularly- 
shaped spots, there are still many uncertainties in the pre- 
cise interpretation of such maps. 

The precise role of, e.g., convection in the outer atmo- 
spheres of red giants is also unclear. Semi-regular variables, 
such as q Orionis, undergo brightness changes of about one- 
half magnitude on a timescale of years; this variability has 
been associated with the intermittent appearance of large 
convection cells on the photosphere. HST observations of 
a Orionis do indeed reveal evidence of a single unresolved 
bright area on the photosphere (Gilliland and Dupree 1996) 
although subsequent spectroscopic observations suggest that 
its origin may be due not to convection but rather to "an 
outwardly propagating shock wave" (Uitenbroek, Dupree 
and Gilliland 1998). 

The study of starspots is clearly, therefore, an important 
aspect of gravitational stellar imaging. Moreover, microlens- 
ing provides a unique opportunity for probing such features 
on the surface of (most probably) slowly rotating red giants, 
which are completely unsuitable candidates for the powerful 
Doppler Imaging technique (see e.g. Strassmeier and Linsky, 
1996, and references therein). 

Finally, we also note that when stars are resolved by 
a lens, the probability of planet detection may be reduced 
since the effect of the finite source is to smear out, and thus 
suppress, the planetary signature (see for instance Vermaak, 
2000, for a recent analysis). Photospheric starspots could 
mimic such features, and hence provide an unexpected noise 
background for planet detection via microlensing, although 
the time of observation of the features and their chromatic 
signature should generally provide a suitable means of dis- 
criminating between spots and planets. 

The structure of this paper is as follows. In §2 we discuss 
the Next Generation stellar atmosphere models used in our 
calculations, in §3 we discuss our model for the microlensing 
of an extended source and in §4 we provide some informa- 
tive examples of the lightcurves produced by spotted stars 
and explore the "detectability" of starspots as a function 
of e.g. size, position and temperature. In §5 we present our 
conclusions. Appendix A describes in detail our geometrical 
treatment of a circular spot. 



2 NEXTGEN MODEL ATMOSPHERES 

Computations of stellar surface brightness profiles have been 
carried out for several decades but until very recently were 
based on an approximate treatment whereby the Planck 
function was used to compute central intensities, 7o, for dif- 
ferent wavebands, and the intensity, I(fi), as a function of 
(the cosine of) emergent angle, fi, was then given by a simple 
parametric model, such as the linear model, namely: 
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7(/x) = / [1-c 1 (1-m)] 



(i) 



The coefficient ci depends on the temperature, gravity and 
chemical composition of the source and on the wavelength 
of observation. Although several improvements have been 
made in recent years (see Valls-Gabaud, 1998, for further 
references) it is clear that some of the underlying hypotheses 
are not realistic. 

The recent 'Next Generation' stellar atmosphere mod- 
els, computed by Hauschildt and collaborators (1999a,b), 
considerably improve upon these simplistic models in several 
important respects. The calculations are carried out assum- 
ing spherical geometry for giant stars, and the intensity pro- 
files are computed directly - without assuming a Planck law 
or a parametric model for the dependence on emergent an- 
gle. Moreover, the intensity calculations are based on a huge 
library of atomic and molecular lines, with about 2 x 10 8 
molecular lines contributing to a typical giant atmosphere 
model at T cff = 3000 K. 

The dramatic difference in the dependence of limb dark- 
ening on emergent angle between the traditional models and 
NextGen models is illustrated in Figure |l|. This figure shows 
the intensity profiles for a giant star of T e ff = 4250K and 
logg = 0.5, in four Johnson colour bands: V, R, I and 
K. The solid curve shows the NextGen profiles, while the 
dashed, dash-dotted and dotted curves denote the linear, 
logarithmic and square root models respectively. We can see 
from this figure that there is a sudden decrease in the in- 
tensity of the NextGen models as one approaches the limb 
of the star - i.e. at [i ~ 0.2. This feature arises from the 
adoption of spherical geometry and the improved modelling 
of molecular scattering in the outer atmosphere of the star, 
and is clearly an effect which one would expect to be highly 
relevant to stellar atmospheres probed by microlensing, but 
is completely absent from the earlier parametric models (e.g. 
Kurucz, 1994) which predict significant intensity all the way 
to n — 0. 

This strong limb feature was detected during the recent 
EROS 2000-BUL-5 event, where a NextGen model was able 
to fit the Ha equivalent width variation across the photo- 
sphere of a K giant (Albrow et al. 2001b), improving the 
earlier calculations by Valls-Gabaud (1996, 1998). Accord- 
ingly the results presented in this paper are calculated us- 
ing NextGen atmosphere models rather than parameterised 
limb darkening laws. 



3 MICROLENSING OF AN EXTENDED 
SPOTTED SOURCE 

The majority of candidate microlensing events are well fitted 
by a point source, point lens model. The characteristic scale 
for microlensing is usually taken to be the Einstein radius 
of the lens, defined as 



Re = 



4GM (Ds-D L )D L 



D S 



(2) 



Here M, Ds, Dl are the lens mass, the distance to 
source and the distance to lens respectively. The angular 
Einstein radius, 6e, is 



0e 



Re 
D L 



4GM (D s - D L ) 



In the case of a point source being lensed by a point 
mass the amplification, A, of the source depends only on 
the projected angular separation, u of the lens and source. 



A(u) = 



m j + 2 
uy/u 2 + 4 



(4) 



and hence is independent of wavelength. 

However when the angular stellar radius is comparable 
to the angular Einstein radius of the lens (i.e. an extended 
source) the integrated, lensed flux from the stellar disk is 
then given by 



F(u) 



I(s,e)A(u)sdsd0 



(•») 



D L D.s 



(3) 



where the source angular radius, R, and angular sep- 
aration, u, are measured in units of 9e- The amplification 
factor, A(u), is simply the point source, point lens amplifi- 
cation given by Equation |^. 

For the background, unspotted photosphere we assume 
that I t (s,8) — I*(s) - i.e. the surface brightness profile is 
radially symmetric, and these values are then taken from the 
NextGen model atmospheres, for the appropriate waveband 
under examination and the effective temperature and logg 
of the source being modelled. 

Our procedure for integrating the flux difference due 
to starspots is very similar to the above, except that the 
appropriate temperature at which we evaluate the surface 
brightness now varies as a function of the polar angle, 6. 
Thus 

AF(u) = [I sp (s,6)-I*(s)]A(u)sdsd9 (6) 

J A 

where A is the projected area of the disk covered by spots. 
We may apply Eq. ^ for the case of several starspots by 
simply adding together the contribution, AF(u)i, from each 
spot - provided only that none of the spots overlaps with any 
other. Furthermore we also apply the NextGen limb dark- 
ening to the spots modelled, i.e. spots close to the limb will 
not be as bright as spots of the same effective temperature 
in the centre of the photosphere. The key problem is, then, 
a purely geometrical one: to determine the (s, 6) coordinates 
marking the boundary of the area covered by each starspot 
at any given time. A similar problem was dealt with by Dor- 
ren (1987), but we describe our method fully in Appendix A. 
The only physical simplification we have made is to assume 
that the effective gravity within each spot is the same as 
the gravity in the photosphere. (It would be computation- 
ally straightforward - but prohibitively time-consuming - 
to allow for variations in effective gravity within each spot). 
Although there is as yet no direct observational evidence to 
support this assumption (Donati, private communication) 
we can show that it is reasonable by the following physical 
argument. Assuming that the growth of spots is regulated 
by horizontal pressure gradients created by the evaporation 
of magnetic flux tubes, Gray (1992) quotes the approximate 
relation, t/s P ot-Ps P ot = constant, between the spot effective 
gravity and pressure - i.e. the larger the spot pressure, the 
lower its effective gravity. Relating P spo t to the stellar mag- 
netic field and inserting typical values for a solar-type star 
with logg = 4.5, one finds log(7 spot ~ 3.6. For giant stars, 
with a smaller magnetic field, the difference in logg should 
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Figure 1. Intensity profiles for a giant star of T e ff = 4250K and logg = 0.5, in four Johnson colour bands: V, R, I and K. The solid 
curve shows the NextGcn profiles, while the dashed, dash-dotted and dotted curves denote the linear, logarithmic and square root models 
respectively. 



be smaller than about 0.2 dex, and so the small effect of 
varying effective gravity within the spot can be neglected. 



4 RESULTS 

4.1 Illustrative examples 

Figures g to [s] present illustrative lightcurves for several 
different configurations of star, lens and spots. In each of 
these figures we have assumed that the angular radius of 
the source is equal to the angular Einstein radius of the 
lens. While this choice of stellar radius is observationally 
unrealistic and produces low amplification events, for illus- 
trative purposes it has the advantage of producing spot fea- 
tures on the lightcurves which are clearly distinguishable by 
eye. We will consider somewhat more realistic values of the 
ratio of stellar radius to Einstein radius later in this sec- 
tion. Lightcurves were computed from the NextGen models 
in the U to K bands; however for clarity only 4 bands are 
shown in the following figures. The V , R, I and K bands 
are represented by a continuous, a dashed, a dotted and 
dashed-dotted line respectively. 

In Figures ^ to [B| three panels are shown. The left hand 
panel shows the stellar disk, illustrating the position and size 
of the starspots and the trajectory (indicated by an arrow) of 
the lens during the event (in fact all four cases are for a lens 
impact parameter, uq — 0). The middle panel indicates the 



(absolute value of the) change, Am, in apparent magnitude 
as a function of time, as given by the formula 



Am 



2.5 log 



( 



-Flsp 
VJusp 



(7) 



where Jlsp and .Ftjsp denote the flux from the lensed, spot- 
ted star and unlensed spotted star respectively. (Note that, 
since the effect of lensing is - of course - to amplify the un- 
lensed flux, the sign of the change in apparent magnitude 
is alwys negative). Similarly the right hand panel shows the 
change, Ams in apparent magnitude resulting specifically 
from the presence of the spot(s), as given by the formula 



Ams = —2.5 log-L 



/-Flsp\ 



(8) 



where -Flsf denotes the lensed flux from the spot-free star. 
(Note that Ams may be either positive or negative). 

While the middle panel indicates the time evolution in 
observed apparent magnitude of the stellar disk plus spot(s) 
during the microlensing event, it is difficult to isolate from 
this panel the contribution to the lightcurve of the spots 
themselves. Thus, although the final panel indicates a flux 
ratio which would not be directly observable, it neverthe- 
less makes clear the residual magnitude change between the 
observed flux from a spotted source and the predicted flux 
from an event with identical stellar and lens parameters (i.e. 
stellar radius, atmospheric limb darkening, lens impact pa- 
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Figure 2. The V, R, I and K microlensing light curves produced by the transit of a lens, with impact parameter uq = 0.06b, across 
the disk of a 5000 K star, of radius 1 0g, with logg = 4.0, with five 4200 K starspots (as shown in the left hand panel) of radius 10°. 
The magnitude changes, Am and Amg are as defined in Equations Q and ^. 
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Figure 3. The V, R, I and K microlensing light curves produced and represented in the same manner as Figure Q but with five hot 
4800 K starspots on the disk of a 4000 K, logg = f.O star. 



rameter etc) but without the non-radial surface brightness 
variations due to the presence of spots. 

Note that Ams as defined by Equation ^ would give a 
wavelength dependent non-zero offset, even in the absence 
of lensing. This arises because the integrated unlensed flux 
from the spotted star would, in any case, introduce a mag- 
nitude change when compared to the unspotted star. Thus 
we can think of Ams as composed of two contributions; 



Amj = AmjL — 2.5 log 




where Ftjsf denotes the unlensed flux from the spot-free 
star. Hence AmsL can be regarded as the magnitude change 
due to the lensing of the spots, while the second term repre- 
sents the non-zero offset discussed above. Although we are, 
of course, primarily interested in AmsL, it is useful to in- 
clude the non-zero offset in the right hand panels of Figures 
^| to ^| as it adds to the clarity of the lightcurve deviations. 

Figure ^| illustrates the V, R, I and K lightcurves pro- 
duced by the transit of a lens across the disk of a 5000 K 
dwarf star and five 4200 K circular spots of radius 10° . These 
spot parameters are motivated by, and are broadly consis- 
tent with, the results of stellar maps derived from Doppler 
imaging studies of active dwarfs (see e.g. Strassmeier and 
Linsky, 1996, and references therein). It can be seen that 
the two cool spots very close to the lens trajectory are the 
only features which contribute significantly to the lightcurve. 
The central spot produces a noticeable 'dip' in the central 



portion of the lightcurve, although the feature becomes less 
prominent at longer wavelengths and is barely detectable in 
the K band. This is not surprising as the contrast in limb 
darkening between 5000 K and 4200 K diminishes at longer 
wavelengths. The spot close to the right hand limb of the 
star produces a signal of slightly smaller amplitude - al- 
though with a similar dependence on wavelength - which is 
due in part to the geometrical foreshortening clearly visible 
in the left hand panel. 

Note also that the position of this spot, in the 'wings' 
of the light curve, renders it very difficult to isolate in the 
middle panel, although its presence is obvious in the right 
hand panel. 

Figure || shows the same lens trajectory and configura- 
tion of five spots as in Figure ^, but now with a temperature 
of 4800 K on a 4000 K giant star. Here the spot radius and 
temperature difference is motivated by the HST observa- 
tions of a Orionis (Gilliland and Dupree 1996). Once again, 
only the two spots which lie very close to the lens trajectory 
are 'imaged' by the lens, producing a peak in the lightcurves, 
which diminishes in amplitude at longer wavelengths. Note 
from the right hand panel that the unlensed offset between 
different wavebands and the peak magnitude change due to 
lensing are both somewhat larger than in Figure ^; i.e. the 
unlensed (and lensed) contrast of hotspots is greater than 
that of cool spots. This is not surprising, since what is im- 
portant is the (log) ratio of spotted to unspotted flux; this 
translates to a larger magnitude change for the hot spot 
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than for the cool spot, although the temperature difference 
has the same absolute value in both examples. 

The lensing signature of multiple spots is also interest- 
ing, particularly if spots are restricted to a given latitudinal 
band. Figures [| and ^ show a 4000 K giant star with a large 
group of 4800 K, 10° radius spots filling a 40° band in lat- 
itude, straddling the equator. In Figure [| the inclination of 
the source star is 90°, so that the lens trajectory intersects 
the spot bands. We can see from the middle and right hand 
panels that the spots produce a complex pattern, with a se- 
ries of 'bumps' each one of which is roughly identifiable with 
an individual spot. The amplitude of each feature, however, 
varies with spot longitude - due to foreshortening and limb 
darkening - and also diminishes considerably with increas- 
ing wavelength as before. 

Another interesting feature is apparent in the right hand 
panel of Figure ^. We can see in the wings of the light curve 
that even in the absence of lensing the presence of the large 
band of spots produces, as one might expect, a substantial 
magnitude offset which is also strongly wavelength depen- 
dent (more than 0.5 mag. at V, reducing to less than 0.2 
mag. at K). In the central part of the lightcurve, as the 
lens transits the star and crosses the band of spots, the flux 
in each waveband is amplified due to lensing. Moreover the 
effect of the lensing is also to increase very slightly the differ- 
ential magnitude offset in the central part of the light curve, 
compared with the offset in the wings. However, this differ- 
ential effect is slightly smaller between the series of 'bumps' 
(when the lens is passing through a 'gap' in the band of 
spots - e.g. at t = to) than it is at the peaks of the bumps 
(when the lens is at its closest to one or more of the spots in 
the band). Thus, the increase in flux due to the presence of 
the spots (which is augmented by the lens but which would 
be apparent even if there were no lensing) is slightly diluted 
when the lens is amplifying the flux from the background 
star more than the spots themselves. 

This interplay between the amplified flux from the spots 
and the background star is more clearly seen in the corre- 
sponding panels of Figure [| which shows the microlensing 
light curves of the same band of spots, but now on a star of 
inclination 60° . Again we see a wavelength-dependent mag- 
nitude offset even when there is no lensing, due to the pres- 
ence of the hot spots. Shortly after the stellar transit begins 
there is a small 'spike' in the light curves, as the lens passes 
very close to one of the spots. Thereafter, however, there 
is a large dip around the central part of each light curve, 
which also slightly reduces the differential magnitude offset 
between the wavebands. This is because, during this time, 
the flux is dominated by the amplified contribution from the 
region of the star very close to the lens, which - as can be 
seen from the left hand panel of Figure a - is spot free. Thus 
we see that even when the lens does not pass directly over a 
spotted region, the presence of spots elsewhere on the star 
still results in a clearly detectable microlensing signature in 
the lightcurve as a whole. 

The systematic deviations from the unspotted light 
curves illustrated in Figures ^ to |E] highlight a potential pit- 
fall of estimating stellar and lens parameters from discretely 
sampled data. In Figure ^, for example, the peak produced 
by the central spot feature could - if sampled with poor 
time resolution - lead to the estimation of a stellar radius 
which was systematically smaller than the true radius. In 



this respect, in addition to the obvious benefit of increasing 
the time resolution of observations, multiwavength photom- 
etry can also be effective in breaking the apparent degener- 
acy between a spotted star and a spot-free star of smaller 
radius. This is because the effect of the starspot is consider- 
ably reduced at longer wavelengths. Similar remarks apply 
to the case of the cool spots on Figure ^, although since 
the amplitude at a given wavelength of the spot feature is 
somewhat smaller for cool spots, the (positive) systematic 
bias in the estimated stellar radius would be less significant 
than the (negative) bias obtained from higher amplification 
hotspots. 

4.2 Investigating spot detectability 

We next carried out an extensive investigation of spot de- 
tectability as a function of spot radius and position and 
lens impact parameter. These calculations were performed 
for the somewhat more realistic case of a star with angu- 
lar radius R, — O.We- This value is comparable to - al- 
beit slightly larger than - the estimated stellar radius for 
MACHO 95-30 (Alcock et al. 1997), and is similar to the 
value adopted in the calculations of HS00. There could, of 
course, be some merit in considering sources of smaller angu- 
lar radius since MACHO 95-30 represents the largest source 
star detected to date. Moreover, a smaller source star would 
also yield a larger microlensing signature. However, given 
the reduced probability of a point lens transit for such a 
star - together with the greatly reduced likelihood of inten- 
sive photometric monitoring during that transit - we adopt 
R, = O.We as a reasonable compromise between consid- 
ering a star with unrealistically large angular radius and a 
star with unrealistically low probability of being observed 
undergoing a point lens transit. (See also Section 5 below). 

Figure [] presents V band contour maps of spot de- 
tectability for a range of spot and lens parameters. In all 
cases, the background stellar temperature and spot temper- 
ature were taken to be 4000 K and 4800 K respectively. 
Rows from top to bottom correspond to lens trajectories 
with impact parameter, uo = 0.0,0.4,0.8 and 1.2, respec- 
tively, normalised to the source radius. Columns from left 
to right correspond to spots of radius 3°, 6°, 9° and 12° re- 
spectively. 

We assume for each light curve that a total of 80 equally 
spaced observations were made during t_E, i.e. the time taken 
for the lens to move 1 angular Einstein radius. (Hence, for 
example t_E = 40 days would imply a sampling rate of two 
observations per day). 

The contours shown denote the contribution to the 
change in magnitude arising specifically from the lensing of 
the spot, defined as AmsL in Equation 9 above. The darkest 
level shown denotes a peak magnitude change of at least 0.1 
mag. (i.e. at least one observation for which AmsL > 0.1); 
the lightest level denotes a peak change of 0.01 mag. 

It is clear from these contour maps that, as one might 
expect, the detectability of starspots is improved greatly in 
regions where the spot lies close to the trajectory of the lens. 

Nevertheless, at least for the larger spots, a substantial 
fraction of the stellar surface would give a detectable signa- 
ture if a starspot were present at that location. For a spot 
of radius 12° and a lens impact parameter uo — 0, for ex- 
ample, we can see from the top right panel of Figure 6 that 
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Figure 4. The V, R, I and microlensing light curves produced by the transit of a lens, with impact parameter uo = 0.0AER, across 
the disk of a 4000 K star, of radius 1 8e, with logg = 1.0, with a band of 4800A* starspots (as shown in the left hand panel) of radius 
10° viewed at an inclination of 90°. The magnitude changes, Am and Am s are as denned in Equations ^ and |^. 
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Figure 5. The V, R, I and K microlensing light curves of the same source as in Figure \\ but now viewed at an inclination of 60° 



the spot would provide a magnitude change of 0.1 mag. or 
greater over about 30% of the visible disk of the star. We 
can compare that with the area covered by the spot itself 
as a fraction of the area of the visible hemisphere, that is 
sin 2 5, so in this case about 4.3%. Hence microlensing is very 
efficient at detecting spots under these conditions. 

It is interesting to note, however, that even for the case 
of uo = 0.129e (i.e. impact parameter slightly greater than 
the stellar radius) one still finds that approximately 40 % of 
the stellar disk yields a magnitude change in excess of 0.01 
mag (compared with the 4% area covered by each spot). 
Although such a signal would clearly be difficult to detect, 
its residuals are comparable to the accuracy levels achieved 
by the best current photometry in microlensing follow up 
surveys. 

The furthest left column of Figure ^| (spot radius of 3°) 
shows, as expected, the smallest 'detectable' regions of the 
photosphere. Indeed, for uq — 0.128e it can be seen that no 
detectable signal would be obtained. Note also the discrete 
features in the top left panel (and in a few other panels) 
which arise because of the small size of the spots and the 
discrete nature of the the temporal sampling. With more 
intensive time coverage these features would become more 
smooth. 

It should be stressed, of course, that the contour maps 
shown in Figure 6 assume that one can predict accurately 
the unspotted lightcurve - i.e. one can accurately esti- 
mate the source radius and lens impact parameter from 
the lightcurves as a whole. Uncertainty in those parame- 



ters would generally impair the sensitivity of the lightcurves 
to the presence of spots by reducing the detectability area 
of the contour maps. A detailed investigation of the impact 
of 'global' parameter errors on the estimation of linear limb 
darkening coefficients from analysis of event OGLE-99-BUL- 
23 was carried out by Albrow et al. (2001a), and a similar 
statistical analysis would be important in the context of 
starspots when confronting real microlensing data. Global 
parameter errors would clearly have the most damaging ef- 
fect for small spots, and indeed might 'wash out' completely 
the signature of 3° spots in the leftmost column of Fig- 
ure 6. In those regions where the spot signature exceeds 
0.1 mag., however, spots would still be easily detectable by 
current folow up surveys, even allowing for uncertainties in 
the global parameters. 



4.3 The feasibility of gravitational imaging 

We have seen in the previous section that high time reso- 
lution observations during a transit event can place useful 
constraints on the existence or otherwise of spots on par- 
ticular regions of the photosphere. The question remains as 
to whether one can constrain the detailed structure of spot 
features from such observations. 

Figure [j] shows a comparison between the light curves 
produced by the event illustrated in Figure [] and a simi- 
lar event, with the same stellar and lens parameters, but 
in which the spots have additional temperature structure - 
specifically a central 'umbra' of temperature 5600 K and ra- 
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Figure 6. The regions of excess magnification of 0.01, 0.05 and 0.1 magnitudes from the lightest gray to darkest due to a spot feature of 
4800 K on a 0.1 logg = 1.0, 4000 K star. Rows from top to bottom correspond to lens trajectories with minimum impact parameter, 
no = 0.0,0.4,0.8 and 1.2, respectively, normalised to the source radius. Columns from left to right correspond to spots of radius 3°, 6°, 9° 
and 12° respectively. 



dius 5° surrounded by a cooler 'penumbra' of temperature 
4800 K (the same temperature as the uniform spots of Fig- 
ure ^|) . The right hand panel of Figure ^ shows the difference, 
in the amplification between the two scenarios, given by 



( Jlsp \ 

no — struc 



(10) 



\ Fusp I struc 

It is evident from Figure Q that the effect of the spot struc- 
ture is small, with a peak deviation from the uniform tem- 
perature case of only AA ~ 0.02, which corresponds to 
Am ~ 0.01 mag. Clearly, then the detection of temperature 
structure, given the accuracy of current photometry, would 
be very difficult. A more serious difficulty, however, is pre- 
sented by the severe ill-posedness of the problem: since spots 
need not in general be circular it is likely that the specific 
photometric signature of temperature structure within a cir- 
cular spot could be closely approximated by a non-circular 
spot of uniform temperature. In essence, one cannot effec- 
tively constrain the 2 dimensional structure of a given spot 
feature from only a 1 dimensional microlensing light curve. 



Similar remarks clearly apply to the photosphere as a whole, 
where groups of (arbitarily shaped) individual spots could 
mimic the signature of a single, larger, spot and vice versa. 

Possible rotation of the source star (which is certainly 
an area of some interest in the study of red giants and su- 
pergiants) is another potential source of degeneracy in the 
signatures of starspots. Depending on the orientation of the 
rotation axis with respect to the lens trajectory, a small 
rotation of the source during the transit of a spot of ra- 
dius, say, 10° might broaden or narrow the spot signature 
in the light curve by an amount equivalent to a change in 
the spot radius of about 1°. Clearly the magnitude of the 
effect depends on the star's rotation period and the lens 
timescale, but the main point is that - without a precise 
model for the shape and size of the spot feature - one can- 
not accurately constrain the rotation period from the spot 
signature. It is interesting to consider whether the use of 
multicolour lightcurves might at least partially break the 
degeneracy between models of spot structure and/or stellar 
rotation. The practical limitation here, however, is the fact 



Microlensing signatures of starspots 9 




-10 1 -2 2 



(t-t )/t E 

Figure 7. A comparison of the V, Ft, I and K microlensing light curves produced by the transit of a lens, with impact parameter 
no = 0.0AER, across the disk of a 4000 K star, of radius 1 AER, with logg = 1.0, with five 4800 K starspots of radius 10° (i.e. Figure 2) 
against an identical source other than additional central temperature structure at 5600 K of 5° (as shown in the left hand panel) within 
the spots. AA is as defined in Equation 10. 



that at longer wavelengths spot signatures are considerably 
supressed and are less sensitive to surface temperature and 
other stellar parameters. In future work, we will, neverthe- 
less, investigate optimal methods for combining multicolour 
data to constrain atmosphere and spot models in this man- 
ner. 



5 DISCUSSION 

The results presented in Section 4 clearly show that the 
impact of non-uniform surface brightness profiles on the mi- 
crolensing light curves of point caustic events is, in princi- 
ple, detectable with current observational precision. There 
are several important caveats, however. 

Firstly, as discussed in Section 4.3, the degenerate na- 
ture of the microlensing light curve, which is a convolution 
over the stellar disk, renders the task of constraining the pre- 
cise number, position, shape, size and temperature contrast 
of starspot(s) an ill-posed problem. Notwithstanding these 
limitations, however, Figure 6 shows that - for a transit or 
near transit event - a significant fraction of the photosphere 
will yield an observable magnitude change due to the pres- 
ence of even a single spot, provided it has sufficient area and 
temperature contrast. (In fact, since the important quantity 
in determining the size of the spot signature is the ratio 
of spotted to unspotted flux, changing the spot tempera- 
ture contrast is a more important factor than changing the 
spot area - at least in the wavebands considered here; for 
the idealised case of black-body radiation and ignoring limb 
darkening, the spot flux will be proportional to the product 
of the area and the fourth power of the temperature). 

Thus it is clear that the failure to detect a significant 
deviation from the light curve signature expected for an 
unspotted source does indeed allow one to place robust lim- 
its on the fraction of the stellar disk (i.e. the 'filling factor') 
which could be covered by starspots. 

Of course the crucial issue here is not simply the pres- 



ence or otherwise of spot features, but also one's ability to 
detect them. 

As we discussed in Section 4.2, an important consider- 
ation here is the question of how well determined are the 
'global' parameters of the event - i.e. the lens minimum im- 
pact parameter and the source radius which determine the 
predicted spot-free lightcurve. These can, in principle, still 
be determined from the overall light curve shape, and in- 
deed the use of multicolour observations - which we have 
shown can be an important diagnostic of the properties of 
the spots - can also improve the accuracy with which the 
global parameters are determined. However, any significant 
uncertainty in the global parameters would weaken the con- 
straints on the spot features since the predicted spot-free 
light curve would be less well determined. A rigorous sta- 
tistical treatment of the impact of errors in the global pa- 
rameters would be a crucial step before robust limits on the 
properties of spots could be inferred from real microlensing 
data. 

The main practical limitation, however, to the use of 
point lens events as probes of stellar surface features is sim- 
ply the trajectory of the lens. Spots which lie in regions of 
the disk far from the lens path will not be readily detected 
under any circumstance. Even regions lying very close to 
the lens trajectory, however, would require excellent pho- 
tometry and high temporal sampling, in order to detect 
smaller features. In the illustrative example of Figure ^| - 
for which we assumed a sampling rate of e.g. approximately 
twice per day for an event with timescale tE — 40 days - 
spots smaller than 3° in radius were undetectable even with 
a minimum detection criterion of having only one observa- 
tion with Am s > 0.1 mag. Clearly a statistically convinc- 
ing detection might require several consecutive data points 
with Am s greater than some chosen threshold. This could be 
achieved either through having larger spots (as is indeed the 
case for the other columns of Figure 6) or by increasing the 
temporal sampling. While the latter is, of course, always de- 
sirable, telescope logistics dictate that - without prior warn- 
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ing that a transit is about to occur - a sampling rate greater 
than about twice per day is not realistic. Transits of a point 
mass lens are rare and even for high amplification events, 
difficult to predict from the pre-transit lightcurve. A more 
'observationally friendly' scenario is the transit of an ex- 
tended source by the caustic structure produced by a binary 
lens. This is because in such a high amplification microlens- 
ing event it becomes necessary to treat every source as an 
extended source and hence the source's surface brightness 
profile can be investigated. Furthermore the intial transit 
'into' the structure acts as a alert, allowing intensive obser- 
vations of the second crossing 'out of the structure to be 
planned. Such 'alert response' photometry is currently being 
carried out with the aim of detecting low-mass companions 
but clearly has scope to image stellar atmospheres, allowing 
the detection or otherwise of photospheric starspots. Gravi- 
tational imaging by non-point lens objects will be addressed 
in a forthcoming paper. 
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APPENDIX A: STARSPOT GEOMETRY AND 
INTEGRATION LIMITS 

Al Coordinate systems 

We calculate the integrated flux from the star, in both the 
lensed and unlensed case, in terms of an integral over the 
projected stellar disk. There are three coordinate systems 
relevant to this calculation: 

(i) (a',cj)') : spherical polar coordinates on the surface of 
the star, with the stellar equator defining <f)' = 0, and with 
a' measured counter-clockwise from the direction which is 
co-planar with the star's rotation axis and the line of sight 
(as shown in Fig. Al). 

(ii) (a, <f>) : spherical polar coordinates on the surface of 
the star, but with polar axis {<f> = -|) defined as the projec- 
tion of the star's rotation axis on the plane of the sky and 
with azimuthal angle, a, measured counter-clockwise from 
the line of sight (as shown in Fig. A2). 

(iii) (s, 6) : projected circular polar coordinates on the 
stellar disk (i.e. the plane of the sky), with 6 measured 
counter-clockwise from the y-axis (see below). 

Figures Al and A2 illustrate these coordinate systems 
and their associated Cartesian coordinate axes. Thus, we 
define the a;-axis to be the line of sight, the z-axis to be the 
projection of the star's rotation axis onto the plane of the 
sky, and the y-axis to be the direction which completes a 
right-handed coordinate system. It is then easy to see that 
the i/-axis and i/'-axis are identical, and the a;'-axis and z'- 
axis are obtained from the z-axis and z-axis by a rotation 
of (-| — i) about the y-axis, where i is the inclination of the 
star. In summary, for a star of radius, R, and a general point 
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Figure Al. Spherical polar coordinates on the surface of the star, 
with the stellar equator defining 0' = 0, and with a' measured 
counter-clockwise from the direction which is co-planar with the 
star's rotation axis and the line of sight. 




Figure A2. Spherical polar coordinates on the surface of the star, 
but now with polar axis {tj> = ^) defined as the projection of the 
star's rotation axis on the plane of the sky and with azimuthal 
angle, a, measured counter-clockwise from the line of sight 



(x, y, z) on the stellar surface 

x = R cos a cos 4> = -R(cos a cos 4>' sin i + sin </>' cos i) 
y = R sin a cos cj> = R sin a' cos </>' 



= Rsm<f) = _R(sin cj>' sin i — cos a cos ( 



cos i 



(Al) 



where all coordinates are expressed in units of the angular 
Einstein radius (AER) of the lens. Consider now the path 
of a point lens, as seen in projection on the sky. Figure A3 
shows the lens trajectory and the position of the lens at 
some general point, L, and at time, t. Here, u denotes the 
impact parameter of the lens and uq denotes the impact 
parameter at the time of closest approach, to, when the lens 
has position angle #o, as shown. Let ts denote the time for 
the lens to move 18 e- The coordinates (j/l, zl), of the lens 
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Figure A3. Lens trajectory 



at position, L, are 

(t - t ) . 
j/l = uo cos do ; sin y 



tE 

Uo sin 0o + — cos 9o 



(A2) 



The projected separation, u, of the lens from an arbitrary 
point, (y,z) — (s cos 8, s sin 9), on the stellar disk is 

u = Uj/l — s cos 9) 2 + (zl — s sin #) 2 ] (A3) 



Consider a general point (x,y,z) on the circumference 
of the spot. We have 

(x — x c ) 2 + (y — y c ) 2 + [z — z c ) 2 = R 2 sin 2 6 



and 

2 . 2 . 2 7~»2 
X + J/ +2 = it 



(A6) 



(A7) 



Combining eqs. A6 and A7 gives, after some manipulation 
R 2 cos 2 8 



xx c + yy c + zz c 



(A8) 



which, as expected, defines a plane perpendicular to the po- 
sition vector (x c ,y c ,z c ). 

A3 Spot visibility conditions 

Consider a unit vector in the direction of the spot centre, 
C. i.e. 

n c — (cos q c cos <j> c , sin a c cos <j> c , sin <j> c ) (A9) 

Let r/ be the angle between the line of sight and n c . Then 
cos 77 = cos ot c cos <f)c- A spot will be fully visible provided 
l]< f -6, i.e. 

cos q c cos <j>c > sin 5 (A10) 
Similarly the spot will be fully invisible provided 
cos a c cos 4>c <— sin 5 (All) 
and partially visible when 

— sin 8 < cos q c cos cj) c < sin 8 (A12) 



A2 Denning the boundary of a circular starspot 

We consider circular starspots - i.e. the locus of points defin- 
ing the boundary of a spot is a small circle of angular radius 
8, say. We suppose that each starspot maintains constant 
radius, temperature and latitude (in the stellar-based coor- 
dinate system) throughout the microlensing event, but that 
its longitude changes if the star is rotating. 

Let (a p ',cj>p') denote the (stellar-based) coordinates of 
the centre of the starspot. If the star is not rotating, these 
coordinates remain fixed; if the star is rotating with period, 
P, and the spot centre transits at time, to, then 

Op'=5(*-*o) ( A4 ) 



We can easily obtain from eqs. Al the projected circular 
polar coordinates, (s p ,9 p ), of P in the observer-based co- 
ordinate system. The circumference of the spot describes a 
planar circle of radius R sin 8. The centre, C, of this circle 
lies inside the star, i.e. 

X c — R COS 8 COS Q c COS (/>c 

y c — R cos 8 sin a c cos </> c = s c cos 9 C 



R cos 8 sin < 



s c sin f 



(A5) 



Note that s c = s p cos 8 and 8 C = 8 P . 

Thus, when seen in projection the centre, C, of the pla- 
nar circle defining the spot boundary is not coincident with 
the centre, P, of the spot on the surface of the star, but does 
lie along the same radial vector joining P to the centre of 
the stellar disk. 



A4 Spot centred on the limb of the star 

Suppose first that x p — 0. It is straightforward to show that 
the spot circumference appears in projection as a straight 
line perpendicular to the radius vector to (y c ,z c ) and the 
integration limits are 

8 c -8 <9 <9 c + 8 <s<R (A13) 



A5 Fully visible spot 

Suppose now that x c 7^ 0. For a fully visible spot, at any 
time the projected spot will appear as an ellipse centred on 
(y c ,z c ). The semi-major axis is perpendicular to the radius 
vector to (y c ,z c ) and some straightforward algebra shows 
that it has length ii = i?sin<5. To determine the semi- minor 
axis we require to solve for the value(s) of s at which the 
spot projection intersects the radius vector through (y c , z c ). 
Clearly, at the points of intersection we have 



y — s cos o c 



z = s sin f 



From eq. hi it follows that 

ZZ C 



R 2 cos 2 8 — yy c 



Combining eqs. A7 and A15 
(R 2 cos 2 8 



yyc 



gives 

22 . 22 



zz c ) + y x c + z x c = K x c 



(A14) 



(A15) 



(A16) 



which, substituting from eqs. A5 and A14, may be reduced 
to the quadratic equation in s 
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R 2 cos 2 S — 2s c s ■ 







(A17) 



This has determinant, A, which some algebra reduces to 
A = 4R 2 cos 2 « c cos 2 4> c sin 8 (A18) 
Hence eq. A17 has roots 

s — s c ± R cos a c cos 4> c sin 5 (A19) 

from which we see immediately that the projected spot el- 
lipse has semi-minor axis h = R cos a c cos <j> c sin 8. 

We can parametrise a general point inside this ellipse 

as 



Ue = uli cos 9e 



ze = UJ h sin 9 e 



(A20) 



Where < u) < 1 and < #e < 27t. The coordinates 
{Ve,Ze) are related to (y, z) via 



y = y E sin ( 
2 = Z£sin( 



+ 2fi cos 6* c + s c cos ( 
— j/_e cos 9 C + s c sin 6 



(A21) 



The integral in eq. A7 may then be expressed in terms 
of u and 9e, viz 



AF = hh 



[I sv {s : e)-h{s)]A{d)ucLjdOE (A22) 



A6 Partially visible spot 

The case where a spot is only partially visible is slightly 
more complicated. Consider the intersection of the projected 
spot ellipse with a circle of radius s on the stellar disk and 
centered on O. Putting y — s cos 8, z — s s in 9, y c — s c cos 9 C , 
z c = s c sin 8 C and substituting in eq. |A8| , gives, after some 
further reduction 



cos(6> - 8 C ) 



R 2 cos 2 8 - XcVR 2 - s 2 



or, writing in terms of a c and 



. ± cos 



R cos 8 — COS Olc COS <f>cV R 2 — s 2 
s(sin 2 etc cos 2 4>c + sin 2 dc) 1 / 2 



(A23) 



(A24) 



Thus, for a partially visible spot we integrate eq. A. 7 
using the limits 



R cos etc cos (h c sin S < s < R 



(A25) 
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